Universal conductivity of graphene in the ultra-relativistic regime 
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We calculate the zero temperature dc conductivity in graphene in the ultimate low-energy regime, 
when retardation effects of the electromagnetic interaction become important and when the full 
Lorentz symmetry emerges. In contrast to what happens with the short range or with the Coulomb 
long-range instantaneous interactions, the dc conductivity is now no longer equal to its non interact- 
ing value, but acquires universal corrections in powers of the fine structure constant. The coefficient 
of the first order correction is computed, and found to be of order one. We also present the result 
for the conductivity in the large-N limit, with A'^ as the number of Dirac fermions species, to the 
order 1/N^. 



PACS numbers: 72.80.Vp, 73.22.Pr,05.10Cc 

INTRODUCTION 

Graphene owes several of its remarkable properties to 
the fact that it admits an effective relativistic quan- 
tum field theory description in terms of massless Dirac 
fermions in tvifo dimension. A dramatic manifestation 
of this fact is in the conductivity properties at half fill- 
ing; recent experiments [1] found that the optical con- 
ductivity in monolayer graphene is essentially constant 
in a wide range of frequencies and very close to the value 
(7r/2)(e^//i) (that is 1/4 in the natural units e = h = 
c — 1 which we will use), which is precisely the value 
found for the system of non-interacting Dirac fermions 
[2]. This result, however, raises a couple of natural ques- 
tions: why the interactions, which are rather strong in 
graphene, do not produce visible corrections to the non 
interacting value of the conductivity? And should the dc 
conductivity in graphene be truly equal to its non inter- 
acting value, or there are in principle corrections which 
may lie inside the experimental errors? 

The computation of the many body interaction effects 
on graphene's conductivity is quite sensitive to regular- 
izations and approximations, and several aspects of it 
have been controversial. In the case of short range in- 
teractions, after first perturbative computations claim- 
ing non vanishing corrections, it was finally rigorously 
proved [3] that there are no interaction corrections to the 
dc conductivity; all the interaction contributions to the 
conductivity cancel out at all orders in the renormalized 
expansion. The exact vanishing of interaction correction 
emerges as a consequence of the Ward identities and the 
irrelevance, in the technical renormalization group (RG) 
sense, of the interaction, [3, 4]. On the other hand, in 
the case of long range Coulomb interactions it has been 
predicted that the dc conductivity is still universal and 



equal to 1/4 [5], the argument this time being based on 
the divergence of the Fermi velocity [6] , and the relative 
(to the kinetic energy) irrelevance of the interaction. The 
low frequency correction to the conductivity was found 
to be not particularly small, due to the slow logarithmic 
increase of the Fermi velocity, which is what is expected 
in absence of (possible) accidental cancellation. A con- 
troversy in the computation of such corrections arose in 
literature, [4, 5, 7-9]. Technically the reason for the con- 
troversy lies in the ambiguities produced by the ultravi- 
olet divergences due to the continuum limit. Recently, 
however, the controversy has been claimed to be settled 
by performing a lattice computation in [10], in favor of 
the value originally found in [4]. 

On the other hand, the Fermi velocity divergence found 
in the Coulomb case at very low frequencies is clearly 
rather unphysical, and simply signals ultimate inade- 
quacy of the usual model of instantaneous Coulomb in- 
teraction. With the increase of the Fermi velocity the 
retardation effects eventually become important, so that 
the retarded current-current interaction must be added 
to the Coulomb density-density interaction. Such effects 
have been analyzed before in [11] , [12] by a RG analysis, 
and it was found that the flow of the Fermi velocity stops 
at the velocity of light c, and, maybe most importantly, 
that the coupling constant (i. e. the charge) in the the- 
ory is exactly marginal. In particular, a lattice model for 
graphene interacting with an electromagnetic field was 
considered in [12], and by iterating the RG it was proved 
that Lorentz symmetry spontaneously emerges, and that 
the system is asymptotically close to the so-called "re- 
duced" quantum electrodynamics (QED^^^). 

The aim of this paper is to compute the corrections to 
the non-interacting value of the universal dc conductivity 
due to the full electromagnetic interaction. We show first 



by exact RG methods that, with the retardation effects 
and with the honeycomb lattice included, in units of e^/ft. 
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where N is the number of four-component Dirac fcrmions 
(A^ = 2 in graphcne), and Ci is an iV- independent con- 
stant. In the ultimate, low-energy regime the conductiv- 
ity is therefore different with respect to its non interact- 
ing vale of 1/4; the dependence on the number of fermion 
components N in the first interaction correction imme- 
diately rules out possible cancellations in general. On 
the other hand, the correction to conductivity, in vac- 
uum with the dielectric constant e = 1 assumed above, is 
still universal, as a consequence of the emerging Lorcntz 
invariance, in the sense that it does not depend on the 
material parameters such as the Fermi velocity, but only 
on the fine structure constant. 

As the lattice model for graphene becomes asymp- 
totically close to QED4,3, which is the fixed point of 
the RG fiow, it is worthwhile computing the dc conduc- 
tivity directly in this continuum model by using stan- 
dard field-theoretical methods. This way an expres- 
sion identical to (1) is found, with the numerical con- 
stant Ci ^ 0.0089319. The leading 0{e'^N) correction 
is identical both in the lattice and in the continuum 
model. In graphene, the number of Dirac fermions is 
N — 2, and this correction is 0.125 and much larger 
than the constant Ci. In any case, as e^ = 47ra, with 
a = 1/137.036 as the fine structure constant, in vac- 
uum we find a/ao = 1 — 0.01064, a small correction, 
and at the moment within the experimental error in [1]. 
The conductivity is therefore in principle, if not in prac- 
tice, different from the non interacting value and, because 
of the smallness of the fine structure constant, with a 
new universal value only slightly reduced from the non- 
interacting one. Finally, we observe in passing that the 
prefactor multiplying the charge e^ in the correction term 
is 0.11607, of the same order in magnitude as found in the 
non-relativistic limit for the static Coulomb interaction 
[4, 5, 10], although of the opposite sign. 



RENORMALIZATION GROUP COMPUTATION 

We can describe graphene by a system of electrons on 
the honeycomb lattice interacting with an electromag- 
netic field. The Hamiltonian is H = H^ + Ha, where 
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with A = [uili + ^2^2), ni = Q, ..,L — \ and ^1.2 = 
i(3,±V3), A = (1,0), h = ^(-l,\/3), 4 - 
2(— 1, — v3) and Ha is the free photon Hamiltonian. 



The physical observables are conveniently obtained in 
terms of the following generating functional: 



where i/'x = ('^x(t'^x+<5 a) ^'^'^ Grassman variables (de- 
noted with a slight abuse of notation with the same sym- 
bol), (T = 1,..,A^, 5j — (0, (5j), P(dtp) is the fermionic 
integration with propagator 



5(x - y) = 
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with X = (xo,x), k = (ko,k), feg = -S-{m + 1/2), k = 
^h + ^62, foi.2 = f (1,±V3), 0^<^ m, < L, |A| = L^, 
V = |i, and il{k) = § Ej=i 2 3 e'^''^'"'*' ^- In the limit 
i ^ 00 1^ Efc ^ S'/e j^, 5 == ^ is the area of the 
exhagonal cell and the integral is over the Brillouin zone. 
The dispersion relation fl{k) vanishes at the two Fermi 



points k^ = (0,f,±^) 



M,^ 



is a Gaussian variable 



while P{da) is the photon integration with propagator 
5f^,iy'w{:x. — y) and w(p) — ^^pf, where x(p) is a cut-off 
functions and the Feynman gauge is assumed. Finally 
the interaction V can be easily deduced from (2) and it 
has the form 



V{a,ip) = '^ e rfx(ao,. 
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where F{a) contains higher order term in a, j^^. 
(.?o,x,.7x)p with 
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j"=l,2,3<T=l 

From the identity WLjiiA, A) = VKl,/3(A + da, Ae'"-) we 
get 

0^-^WLAA + da,Xe'''-) = 0, (7) 

and the derivatives of the above relations provide 
a set of Ward identities. Calling iir^^y(p) = 

4 --, ^;fi ' ' |o, p — (w,p), the dc conductivity at 

zero frequency is defined via Kubo formula a = 
\imui^o~e^—Ki^i{uj,0). Note that Kij{p) is the sum 
of the truncated current-current correlation and of the 
diamagnetic term, which is a constant in p; by the 
Ward identity obtained from (7) with a derivative with 
respect to A, we get limp^^o lim^^o -ft'ii(p)|p2=o = 



limpj_>.olim^_).o ^-?iroi(p) = and, as Kii{p) is contin- 
uous at weak coupling, we can reverse the limits so that 
the conductivity can be written as 



a = lim 



,K,,ioj,0) - Ku{0,0) 



(8) 



There is therefore no need of computing the diamagnetic 
term, but it is sufficient to compute the current-current 
correlation and subtract the value in p = (0,0). 



see [12]. The starting point consists of writing the Grass- 
man variables as sums of variables with momenta closer 



and closer to the two Fermi points k^ = (0, ^,±:^y=), 



that is 
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where tpe lives on a shell of momenta distant 0(2'') 



from k| 



Similarly, we can write a^ = ^^^^ „ a'^ 



h' 



and 



a^ has propagator //i(k)w(k) with //i(k) non vanishing 



Kij{p) can be computed by Wilsonian RG; for details, for 2 < |k 



(^i,ai),..(^ 

tional can be written as 



< 2'*+^. After the integration of the fields 
a'*), one finds that the generating func- 



J 



^Wl,;i(A,\) ^ B''(A,\) 



I n P(rf^i-''^) f Pida^^'')e''"^^ 



e=±l 



\p(<'-),a'-^'-)+AA) 



(10) 



where P{d'4>l ) has the propagator 



ffl'^^k' 






Vh^{k) iko 



(11) 



Zh is the wave function renormalization and Vh is the ef- 
fective Fermi velocity. Note also that (12) can be written 
as 



Jh) 
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so that it becomes ever closer to the Dirac propagator. 
By power counting, the scaling dimension of the inter- 
actions is D = 3 — n if n is the number of fields; the 
local terms quadratic in a, describing the photon mass, 
are absent by exploiting the Ward identity generated by 
(7) (see App. El of [12] for an explicit computation at 
one loop ) while the marginal terms a^da are vanishing 
by symmetry; in the same way the local terms ■0"'"^ are 
vanishing by parity, and ip~^dip contribute to the wave 
function renormalization and the Fermi velocity. The ef- 
fective potential is given by (at A = for definiteness): 



/i X 



where Fh are the irrelevant terms, with negative dimen- 
sion D < 0. 



From the Ward identities (7) we get 



Z, 



O.h 



Zh 



= l + 0(e2). 



ZhVh 



= l + 0{e^). (14) 



Moreover, by symmetry Zi^h = Z — 2,H and by an ex- 
plicit computation, 

3^2-"'' ,^;^^l + A(l-^;)2^^ (15) 

with ri,T] > and O(e^) and A is a bounded function. 
Therefore the wave function renormalization is diverging 
(i. e. there is an anomalous dimension), and by iterating 
the RG the Fermi velocity converges to the velocity of 
light. 

The model considered in [11] therefore emerges nat- 
urally starting from the lattice model. Moreover, the 
Lorentz symmetry is restored in the infrared limit, as 
the Fermi velocity flows up to the velocity of light (which 
in our units is 1). The effective couplings e^i^o — s~§^ 
^•'^ flow to a line of fixed points 
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(16) 



that is, the theory is exactly marginal. One obtains then 
a renormalized expansion for the current-current corre- 
lations in terms of the effective couplings e^^t, namely 
(7 = (t(°) + (T^^) + .... This differs from perturbative ex- 
pansion in that there are no infrared divergences and all 
coefficients are bounded. We get 
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2tt Je (2^)2 ZhZh 



.Tr{aigi'')(k')ai [gW(k' + (c, 0)) - g('')(k')] }, (17) 



where aj are Pauli matrices and iV = 2 in graphene. 
Note the presence of a factor -^ for any fermionic hne, 
and of a factor Zi^ for any vertex. The presence of such 
factors could radically alter the conductivity properties; 
for instance if we do not take into account the vertex 
renormalization, that is we replace Zij-^ with unity, one 
would get CT*^°^ = 0. On the contrary, thanks to the Ward 



identity and (14) we can replace 



(z. 



ZhZh 



with Vh up to 



0{e ) terms. The integral then still appears to be non 
universal, as it depends of the effective Fermi velocity 
and on the lattice details. However it is not so; we can 
write integral over the Brillouin zone B as sum of two 
terms, one |rj(fc)| > £ and the other |ri(A:)| > e; the for- 
mer is uniformly convergent as cj — > 0^: therefore, we can 



J 



exchange the integral with the limit and check that the 
integral of the limit is zero simply because the integrand 
is odd in Uq. In the latter we use (12) neglecting the cor- 
rections (as the size of the integral is arbitrarily small) ; 
the dependence from Vh disappears through a change of 
variables and we finally get a'^'^'^ — ^^T +0(iVe^). More- 
over, as shown in [13], the low frequency corrections are 
0(e2w2)^ 

Let us consider now a'^'^^] there are three possible con- 
tributions, cr'^^ ^ (ya +(y\, + Cc but only one of them, 
which we call Uc , is proportional to N"^. Therefore, if 
we find it to be non vanishing we can safely conclude that 
the dc conductivity is different from the non-interacting 
value, at least for a generic N . The value of such term is 
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where h{u}) is such that 2'^'"^ ^ |u;|. In writing the above 
expression we have used that only the renormalized parts 
of the to "bubble" diagrams contribute, as their local part 
is vanishing. Again using that Zi^h/Zh — v^ + O(e^) we 
get ai^'' ^ -(eV2)(iV/8)2 + 0{e^N), and finally (1) is 
found. The A^-dependence allows us to exclude the pos- 
sibility of cancellations, and we can conclude that the dc 
conductivity is different from its non interacting value, 
with a leading correction which is universal. Note the 
difference with the case or Coulomb interactions: in such 
a case at one loop Vh ^ VI^L while Z^, Zi,h,ei,h are 

(2) 

essentially constants so that CTc vanishes. Similarly, for 
Hubbard interactions the photon propagator l/2u! should 



be replaced by a constant, and again a, 
ishing. 



(2) 



would be van- 



EFFECTIVE QED DESCRIPTION 

The previous analysis shows that lattice graphene sys- 
tem flows, by iterating the RG, to a fixed point expressed 
by the the continuum QED^^j,. It is useful then to study 
the conductivity properties in the continuum theory (as- 
suming as usual that the value of the conductivity de- 
pends only on the fixed point of the RG,[14]) in order to 
get more information on the subleading corrections. The 
action for the effective model in 2 -I- 1 dimensions is 



S = j dx[*,(x)7^(9^ - e(a^(x) - yl4x)))*,(x) 
+ - / d-xdyW-^^ - y)a^(x)a,(y)], (19) 



where ^ is a four-component fermionic field, i = 1, ...N , 
x = {xQ,Xl,X2)^ and the summation convention is as- 
sumed. Note that x = (xi,X2) is here a continuum vari- 
able while is the previous section was a site on the hon- 
eycomb lattice. The bare gauge field propagator is 
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with the usual transverse projector 

n^i>(q) = s^^ — ^ 



(H^.(q)+/3^), (20) 
q^ 



(21) 



The above field theory is closely related to the three 
dimensional quantum electrodynamics {QED3),[15] with 
one important caveat: the "Maxwell" term is now al- 
ready at the bare level non-analytic in momentum, and 
proportional to |q|, and not to the usual q2. It can 
be obtained from the reduced quantum electrodynam- 
ics QED42 in which the electromagnetic fields live in 
3-1-1 dimensions, but are coupled to fermions which are 
confined to the lower, 2 + 1-dimensional "brane" , by "in- 
tegrating out" the out-of-plane components of the vector 
potential. [16] This procedure in general also changes the 
effective gauge-fixing parameter from /3' in the original 
3 -|- 1-dimensional theory into /3 in the above expression, 
as in /? = {l+/3')/2. We see that only the Feynman gauge 
(/? = 1) remains invariant under this dimensional reduc- 
tion, which makes it the most convenient one from the 
practical point of view. The non-analyticity of the gauge 
field propagator around q = can be understood as the 
reason for the exact marginality of the charge coupling. 
[17] 



Fermions appear quadratically in the action and can 
be (formally) integrated out. If we redefine the fields as 
ea^ — > a^, and eA^ — > A^, and then shift the fluctuating 



field as Qfi — Af^ — > a^, the result of this integration would 
be the action 
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The first term proportional to N is the familiar one-loop 
polarization in the QED^\19\ and the second term is the 
quasi-Maxwell term, which now after the shift of vari- 
ables also includes the external probe; V{a) is a sum 
of monomials in the fluctuating gauge field a with de- 
gree > 4. Note that V{a) docs not contain the external 



probe A , which after the shift appears only in the quasi- 
Maxwell term. This allows one to differentiate with re- 
spect to the external probe, and so to obtain the current- 
current correlation function in terms of the exact gauge 
field propagator, DJ^" = (a^ai/): 
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Current conservation, on the other hand, dictates that 
the exact gauge-field propagator has the form: [20] 



Z)-(q) = i^(An,.(q)+/3^5^; 



(24) 



where the dc conductivity a, in units of e^/ft, is simply 



^=1(1_M). 



(26) 



where ^ is a function of the number of fermions N and of 

the coupling e^. Inserting this form into the expression Note that the current-current correlation function is com- 



for the current-current correlation function, we find that 

OM(q)j-(-q)> = ^|q|nM.(q), (25) 



pletely independent of the gauge-fixing parameter /3, just 
as one expects. We can rewrite (22) as a perturbed Gaus- 
sian action 
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||q|(n^.(q) + iS^)((2a^(q)A.(-q) + A4q)A.(-q))}} + Via). 
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where 7 = {2/(3e^){l/{N/8 + 2/e'^). In the gaussian ap- 
proximation, that is neglecting the higher-order terms 
given by V{a), the functional integral can be explic- 
itly performed and the constant A is given by Aq — 
{{N/8) + (2/e2))-i. From (26) we this way find the loffe- 
Larkin-like [21] result for the conductivity 



which can be interpreted as the addition of the fermion's 
and the gauge-field's resistivities into the total resistivity. 
Expanding to the first power in the weak charge coupling 
e^ yields a = {N/8){l — {Ne^/16) + ..), in agreement with 
the lattice computation. 
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To go beyond the gaussian approximation we can ex- 
v2°j pand in powers of the effective charge in the theory 



eV(2 + (A^eV8)); one finds 
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where x = (92 - 9T:^)/{{4Trfl8) [18]. In the weak cou- 
phng regime e^ <C 1/N, after expanding in powers of 
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(30) 



so that Eq. (1) is recovered, with the numerical value of 
the constant 

23 1 

Ci = 8x = — — - - - = 0.0089319. (31) 

(Jtt)^ 4 

On the other hand in the large-N limit, N 3> 1/e^, we 
expand in powers of 1/iV, 



2 , 2 
'^=-^(1 



1 



e2f + ^+8x(l 



16 



) + 



(32) 



so that 



- = |(l-]^ + ]S(2 + Ci^^) + 0(A^-^))- (33) 



CONCLUSION 

In conclusion, we have presented the computation of 
the universal zero temperature dc conductivity in the ul- 
timate, relativistic, low-energy regime in graphene, when 
the Fermi velocity has reached the velocity of light. Al- 
though this regime lies beyond presently available ex- 
perimental conditions, such as the temperature and the 
sample sizes, the issue of ultimate value of the dc con- 
ductivity is theoretically interesting, and presents an im- 
portant question of principle. We find that the ultimate 
value of the conductivity is universal, and dependent only 
on the fine-structure constant of the media surrounding 
the graphene sheet, but in principle different from the 
non-interacting value seen in the experiment. The differ- 
ence from the non interacting value is of the order of the 
fine structure constant itself, however, and therefore of 
the relative size of the order of one percent. 
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